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We present results for the nucleon electromagnetic form factors, including the momentum transfer 
dependence and derived quantities (charge radii and magnetic moment). The analysis is performed 
using 0{a) improved Wilson fermions in Nf = 2 QCD measured on the CLS ensembles. Particular 
focus is placed on a systematic evaluation of the influence of excited states in three-point correlation 
functions, which lead to a biased evaluation, if not accounted for correctly. We argue that the use of 
summed operator insertions and fit ansatze including excited states allow us to suppress and control 
this effect. We employ a novel method to perform joint chiral and continuum extrapolations, by 
fitting the form factors directly to the expressions of covariant baryonic chiral effective field theory. 
The final results for the charge radii and magnetic moment from our lattice calculations include, 
for the first time, a full error budget. We find that our estimates are compatible with experimental 
results within their overall uncertainties. 
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I. INTRODUCTION 

The electromagnetic form factors, Ge and Gm, of the 
nucleon encode information on the distribution of charge 
and magnetization and are among the key quantities de¬ 
scribing its internal structure. Experimental measure¬ 
ments of these quantities in ep scattering processes have 
a long history (see, e.g. the review [T]) and have been 
pushed to ever higher precision [1 In spite of the 
fact that nucleon electromagnetic form factors have been 
studied extensively in theory and experiment, there are 
several open questions. The first concerns the deviation 
between the ratio Ge/Gm as determined using the tra¬ 
ditional Rosenbluth separation technique and the result 
obtained from recoil polarization gHll at squared mo¬ 
mentum transfers larger than 1 GeV^. Secondly, 
prompted by the observed discrepancy between the pro¬ 
ton charge radius extracted from the Lamb shift in 
muonic hydrogen mis] and the value obtained by using 
the electron as a probe [31 [ID], there is a strong interest 
in new experimental measurements of form factors in the 
regime of very small Q^, as well as in further theoretical 
studies, in order to reduce the inherent systematics. The 
third open issue concerns our understanding of the inter¬ 
nal structure of the nucleon in terms of the underlying 
gauge theory of QGD. Nucleon form factors have been 
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studied extensively in simulations of QCD on a space- 
time lattice [IIH23I, and although these calculations are 
quite straightforward, they mostly fail in reproducing the 
experimentally observed Q^-dependence of Ge and Gm- 
As a consequence, lattice estimates for the electric charge 
radius derived from the slope of Ge at vanishing are 
typically underestimated compared to the results derived 
from ep scattering data. It is widely believed that system¬ 
atic errors in lattice calculations must be held responsible 
for this deviation. 

In addition to systematic errors induced by non-zero 
lattice spacings, finite volumes and by uncertainties asso¬ 
ciated with the chiral extrapolation, the issue of contam¬ 
ination from excited states in calculations of nucleon cor¬ 
relation functions has recently come to the fore as a pos¬ 
sible explanation for the deviation between experimental 
and lattice estimates of the electric charge radius. 

In this paper we present a detailed investigation of sys¬ 
tematic effects in lattice calculations of nucleon form fac¬ 
tors arising from excited state contributions. In partic¬ 
ular, we apply the technique of summed operator inser¬ 
tions |D3H3Sj which has proved very useful in our earlier 
calculation of the axial charge of the nucleon |DZ]. Fur¬ 
thermore, we address in detail the chiral extrapolation to 
the physical pion mass, by employing several variants of 
baryonic Chiral Perturbation Theory (ChPT). 

Our simulations are performed in two-flavour QCD 
with a mass-degenerate doublet of up- and down-quarks. 
Since excited state contamination is an issue for lattice 
simulations with any number of dynamical quarks, the 
question whether estimates for nucleon charge radii and 
magnetic moments may be biased can be adequately ad¬ 
dressed in this set-up. There is ample evidence [28] that 
there are no discernible differences between QCD with 
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N{ = 2 and = 2 + 1 flavours at the few-percent level. 
Therefore, the observed deviation between lattice QCD 
and experiment is far too large to be explained by the 
presence or absence of a dynamical strange quark. 

Our central findings include the observation that 
excited-state contaminations have a sizeable influence 
on the form factors extracted from the still widely-used 
plateau method applied to ratios of three- and two- 
point functions at least up to source-sink separations of 
^ 1.5 fm. The use of summed insertions, while gener¬ 
ally an important tool in suppressing excited-state effects 
on hadron structure quantities, cannot reliably exclude a 
residual bias, in particular when comparing with the re¬ 
sults of fits which include excited states explicitly. More¬ 
over, for the first time, we apply the full framework of 
covariant baryonic chiral perturbation theory to 

the simultaneous determination of the form factors near 
= 0 and at the physical pion mass. From a careful 
study of all relevant systematic effects, we are able to give 
a full error budget. Our final results for various charge 
radii and the anomalous magnetic moments k are listed 
in eq. (341 below. We observe agreement with experiment 
within the accuracy of our calculation, including system¬ 
atic errors. However, the overall uncertainty is too large 
to have an impact on the proton radius puzzle. 

This paper is organized as follows: In section|nj we de¬ 
scribe our lattice setup, including details of the ensembles 
used and observables measured, as well as our evaluation 
of statistical errors. In section 11111 we discuss the analy¬ 
sis methods we employed to study and suppress excited- 
state contributions. The Q^-dependence of the measured 
form factors, and the values of the charge radii and mag¬ 
netic moment determined from dipole fits on each ensem¬ 
ble, are presented in section |IV] In section [Vj we discuss 
in detail the chiral fits to the form factors which we use 
to obtain our final results. Section Ivll contains our con¬ 
clusions and a brief outlook. 

A discussion of the impact of the use of Lorentz non¬ 
covariant interpolating operators obtained from smear¬ 
ing the quark fields in the spatial directions only on the 
Lorentz invariance of the results so obtained is contained 
in appendix For ease of reference, we provide tables 
containing the full set of our results for the form factors 
at all values of on all ensembles in appendix [b| 


II. LATTICE SETUP 


reads 


{Nip',s')\V,^mN{p,s)) = 


u{p', s') 




I - 

2mN 


F2{q^) 


u[p, s), 


( 2 ) 


where u{p, s) is a Dirac spinor with spin s and momentum 
p, is a Dirac matrix, | [ 7 '', 7 "], and ton denotes 

the nucleon mass. 

The four-momentum transfer q = p' — p is expressed 
in terms of the energies and three-momenta of the initial 
and final states as 

q^ = -Q^ = ^Ep,-Epf-{p'-pf. (3) 

In this paper, we focus on the iso-vector form factors. 
By assuming isospin symmetry, one can show via a sim¬ 
ple application of the Wigner-Eckart theorem applied in 
isospin space that 

(p(p',s')|u7^u-J7^d|p(p,s)) = 

(4) 

(p(p', s') 114 '^ |p(p, s)) - (n(p', s') 1144 |n(p, s)) , 

where |p) and |n) refer to one-proton and one-neutron 
states, respectively. The expression on the left-hand side 
is suitable for lattice QCD calculations, while the right- 
hand side allows one to compare the results to experi¬ 
mental measurements. 

The Dirac and Pauli form factors give rise to the 
helicity-preserving and helicity-flipping contributions to 
the amplitude, respectively. The electric and magnetic 
(Sachs) form factors Ge and Gm are obtained as linear 
combinations of Fi and F 2 , 

GE(g") = i"i(g") + ^F2(g"), (5) 

GM(g') = Fi(g2) + F2(g2). (6) 

They can be determined from ep scattering experiments 
by decomposing the measured differential cross section 
through the Rosenbluth formula [33]. The form fac¬ 
tors may be Taylor-expanded in the squared momentum 
transfer g^, 

GE,M(g^) = Ge,m( 0) -f g(?'E,M)9^ + G(g‘^)^ , (7) 


A. Observables and correlators 


from which the charge radii of the nucleon may be deter¬ 
mined: 


The matrix element of the electromagnetic current 

yera=\ui^U-^dYd+ ... ( 1 ) 

between one-nucleon states can be expressed in terms 
of the Dirac and Pauli form factors, Fi and F 2 . In 
Minkowski space notation the form factor decomposition 


'E,M 


) = 


6 9GE,M(g^) 


GE.M(g^) i9g2 


( 8 ) 




Electric charge conservation implies Ge( 0) = 1, while 
the magnetic moment p of the nucleon, in units of the 
nuclear magneton e/2TON, is obtained from the magnetic 
form factor at vanishing g^, Gm(0) = p. 
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FIG. 1. Baryonic two-point and three-point functions (left 
and right panels respectively). 


From the asymptotic behaviour of one can 

then extract Ge and Gm for space-like momenta = 
—q^ > 0 via 


Rvo{q,t,ts) 


TTljSf Eq 


2En 


Gr™(Q"), (14) 


and 


Lattice simulations allow for the determination of 
hadronic matrix elements by computing Euclidean cor¬ 
relation functions of local composite operators]^ To this 
end, one considers the nucleon two-point function 


Rei?y,(q,t,G) 




2Eq{Eq + tun) 




(15) 


C2{p,t) = (9) 

X 

where fE'“(a;, t) denotes a standard interpolating operator 
for the nucleon, and F is a projection matrix in spinor 
space. The three-point function of a generic (Euclidean) 
vector current is given by 


where the e-symbol in the last equation denotes the an¬ 
tisymmetric tensor with ei 23 = +1, and the superscripts 
“bare” remind us that, in general, the vector current re¬ 
quires renormalization in the lattice-regularized theory. 

B. Simulation details 


G3,y,(q,<,G) = ^e*'^'^r^„(vI/“(a;,G)F^(y,t)vI/''(0)), 

x,y 

( 10 ) 

where q = p' — P- For the nucleon correlation functions 
considered in this work, the projection matrix 

r = ^(1+ 7o)(i+ * 7573 ) (11) 

ensures the correct parity of the created states and gives 
the nucleon a polarization in the z-direction, which is re¬ 
quired to extract the magnetic form factor. In the above 
expression for the three-point function the vector current 
is inserted at Euclidean time t, while the Euclidean time 
separation between the initial and final nucleons is de¬ 
noted by tg. Eigurej^ shows the corresponding diagrams 
of the two- and three-point functions. Note that for the 
iso-vector vector current considered in this work, quark- 
disconnected diagrams cancel. Moreover, our kinematics 
is chosen such that the final nucleon is always at rest, i.e. 

p' = 0, q = -p. (12) 

The electric and magnetic form factors are easily deter¬ 
mined from suitable ratios of correlation functions. Here 
we follow ref. [34] and use the ratio found to be most 
effective in isolating the desired matrix element. For our 
chosen kinematics it reads 

Rv^iq,t,tg) = 

C3y^{q,t,tg) l c2{q,tg-t)C2{0,t)C2{0,tg) (13) 

G2(0,G) y 02 ( 0 , ts-t)C2{q,t)C2iq,ts) ' 


^ From here on we use Euclidean notation for the position- and 
momentum-space vectors, as well as for the Dirac matrices. 


Our calculations have been performed on a set of en¬ 
sembles with Nf — 2 flavours of G(a)-improved Wilson 
quarks and the Wilson plaquette action. For the improve¬ 
ment coefficient Csw we used the non-perturbative deter¬ 
mination of ref. |35] . The gauge configurations have been 
generated as part of the CLS (Coordinated Lattice Sim¬ 
ulations) initiative, using the deflation-accelerated DD- 
HMC [^ [37] and MP-HMC [38] algorithms. Table |!] 
provides details of the lattice ensembles used. 

For the calculation of three-point correlation functions 
we employed the point-split iso-vector current 

= I + A)(i+ 

- V'(a;)(l - -i^)Uf,{x)T^ip{x + /})) , 
as well as the local vector current 

(17) 

Here, tp denotes an isospin doublet of up and down quark 
fields, and is the Pauli matrix acting in isospin space. 
While the point-split current is conserved and satisfies 
the corresponding Ward identity, the local vector current 
must be renormalized. The expression for the renormal¬ 
ized current in the 0{a) improved theory reads [33] 

= Zv(l + 6va™q)(Pj°“ + acyd.T ^,), (18) 

where Wq denotes the bare subtracted quark mass, by 
and cy are improvement coefficients, and T^''{x) = 
—'0(x)4[7^,7y]r^'0(x) is the tensor density. We used the 
non-perturbative estimate for the renormalization factor 
Zy in the two-flavour theory of refs. [dO] [41]. On the 
other hand, the conserved vector current, while not sub¬ 
ject to renormalization, requires 0{a) improvement even 
at tree level. In this work we neither used the improved 
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Run 

L/a 

d 

n 

am-n- 

arriN 

rtiTrL 

Nctg 

^meas 

a [fm] 

m-n [MeV] 

ra-K/mN 

A3 

32 

5.20 

0.13580 

0.1893(6) 

0.546(7) 

6.0 

133 

2128 

0.079 

473 

0.346(5) 

A4 

32 

5.20 

0.13590 

0.1459(7) 

0.488(13) 

4.7 

200 

3200 

0.079 

364 

0.299(7) 

A5 

32 

5.20 

0.13594 

0.1265(8) 

0.468(7) 

4.0 

250 

4000 

0.079 

316 

0.270(5) 

B6 

48 

5.20 

0.13597 

0.1073(7) 

0.444(5) 

5.0 

159 

2544 

0.079 

268 

0.242(3) 

E5 

32 

5.30 

0.13625 

0.1458(3) 

0.441(4) 

4.7 

1000 

4000 

0.063 

457 

0.330(3) 

F6 

48 

5.30 

0.13635 

0.1036(3) 

0.382(4) 

5.0 

300 

3600 

0.063 

324 

0.271(3) 

F7 

48 

5.30 

0.13638 

0.0885(3) 

0.367(5) 

4.2 

250 

3000 

0.063 

277 

0.241(4) 

G8 

64 

5.30 

0.13642 

0.0617(3) 

0.352(6) 

4.0 

348 

4176 

0.063 

193 

0.175(3) 

N5 

48 

5.50 

0.13660 

0.1086(2) 

0.329(2) 

5.2 

477 

1908 

0.050 

429 

0.330(2) 

N6 

48 

5.50 

0.13667 

0.0838(2) 

0.297(3) 

4.0 

946 

3784 

0.050 

331 

0.283(3) 

07 

64 

5.50 

0.13671 

0.0660(1) 

0.271(4) 

4.4 

490 

1960 

0.050 

261 

0.244(3) 


TABLE I. Details of the lattice ensembles used in this study, showing the lattice extent, L, where T = 2L; the values of the 
bare parameters /3 and k in the lattice action; the pion and nucleon masses (am,r and aniN)', the number of measurements, 
Nmeas = A^cfg X Nsrc', the lattice spacing, a; the pion mass, m,r, in physical units, and the ratio mnlmN of the pion and nucleon 
masses. 


version of the point-split vector current, nor did we com¬ 
pute matrix elements containing the derivative of the ten¬ 
sor current. Therefore, our results for form factors and 
charge radii are not fully 0{a) improved; hence, neglect¬ 
ing the by term in eq. (18) is consistent. 

The interpolating field for the proton was chosen as 


malism intact, breaks the relativistic covariance of the in¬ 
terpolating fields constructed from smeared quarks. This 
issue has not been studied previously in any great detail 
in the context of nucleon form factors. In appendix [X| we 
give a brief explanation why the relativistic invariance of 
our results is not affected. 


4'“(a;) = €abc {ul{x)C^ 5 db{x)) u“(a;), (19) 

with Gaussian-smeared quark fields |32] 

'^ = {1 + kgA)^ Ip , ( 20 ) 

where the links in the three-dimensional covariant Lapla- 
cian A were APE-smeared [43] in the spatial directions 
to further enhance the projection properties onto the 
ground state and help reduce the gauge noise. Correla¬ 
tion functions were constructed using identically smeared 
interpolating fields at both source and sink to ensure that 
the two-point functions are given by a sum of exponen¬ 
tials with positive coefficients. The smearing pa¬ 

rameter KG and the iteration number N were tuned so 
as to maximize the length of the effective mass plateaux 
in a variety of channels. A widely used measure for the 
spatial extent of a smeared source vector is the “smear¬ 
ing radius” rgm (for a definition see e.g. eq. (2.6) in |44]f. 
We note that our choice of kg and N corresponds to 
Tsm ~ 0.5 fm. As was first noted in [33|, the standard 
Gaussian smearing procedure becomes rapidly ineffective 
for baryons as the lattice spacing is decreased. Alterna¬ 
tively one may employ “free-form smearing” |44j which, 
however, cannot be readily applied at the sink. There¬ 
fore, all results presented in this paper have been ob¬ 
tained using standard Gaussian smearing at both the 
source and sink. Note that we did not employ boosted 
Gaussian smearing [3S| either, because the boost is small 
for the nucleon, and the gain in terms of an enhanced pro¬ 
jection on the ground state is expected to be marginal. 

Smearing the quark fields in the spatial directions only, 
while required in order to keep the transfer matrix for¬ 


To compute the three-point function, we use the “fixed- 
sink” method, which requires an additional inversion for 
each value of ts, but allows both the operator inser¬ 
tion and the momentum transfer to be varied without 
additional inversions [46] . In order to realize a range 
of values for the squared four-momentum transfer 
we have computed the two- and three-point correlation 
functions for several spatial momenta q = n27r/L, with 
|np = 0,1,2,...,6. 


The ratios Ry^ of eq. (131 contain a particular combi¬ 
nation of nucleon two-point functions, in order to isolate 
the relevant matrix element. In our analysis the two- 
point functions which enter i?vy were represented by sin¬ 
gle exponential fits. For non-vanishing momenta q, the 
nucleon energies were determined from the nucleon mass 
using the continuum dispersion relation. Compared to 
determining the nucleon energies directly from the ex¬ 
ponential fall-off of C 2 {q,t), we found that this proce¬ 
dure resulted in smaller statistical errors, whilst produc¬ 
ing compatible results. 


In order to express dimensionful quantities in physical 
units, we determined the lattice spacing for all our ensem¬ 
bles using the mass of the If baryon, as described in m- 
More recently, the ALPHA collaboration has published 
accurate values for the lattice spacing determined from 
the kaon decay constant, Jk [IHIIISI- While the central 
values differ slightly, both determinations are well com¬ 
patible within the quoted uncertainties. We have verified 
that uncertainties in the scale-setting procedure have no 
significant influence on the values of the charge radii in 
physical units. 
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C. Statistics and error analysis 

We computed two- and three-point correlation func¬ 
tions on all ensembles listed in Table HI In order to in¬ 
crease statistics, we used multiple sources spread evenly 
across the lattice on each gauge configuration. The total 
number of measurements for each ensemble is listed in 
Table |l] Statistical errors were estimated using a boot¬ 
strap procedure with 10,000 bootstrap samples. 

Simulations at the fine lattice spacings considered here 
are known to be affected by the critical slowing-down of 
the smooth modes of the gauge field and the freezing of 
the topological charge [50], leading to potentially long 
autocorrelation times. Ignoring the long tails in the au¬ 
tocorrelation function may lead to a significant underes¬ 
timation of statistical errors m- Since the correlation 
functions of the nucleon studied here are intrinsically very 
noisy, however, one may expect that the contributions 
from the tails have relatively little influence on the over¬ 
all statistical error. 

We have investigated the impact of autocorrelations 
on our results by performing a binning analysis prior to 
applying the bootstrap procedure. To this end we fo¬ 
cussed on the N6 ensemble, which is based on a long 
Monte Carlo sequence, comprising 8040 molecular dy¬ 
namics units (MDUs) in total. Our findings indicate only 
a marginal increase in the statistical error of the electric 
form factor and the nucleon mass, which amounts to 2% 
at most. We conclude that, for the purpose of comput¬ 
ing nucleon hadronic matrix elements and masses, our 
ensembles are sufficiently decorrelated. 


III. EXCITED-STATE SYSTEMATICS 


“bare” on the form factors, since we assume that the 
ratios have been appropriately renormalized. 

It is well known that nucleon correlation functions suf¬ 
fer from an exponentially increasing noise-to-signal ratio 
[521153] , which imposes a limit on the source-sink separa¬ 
tion ts which can be realized with reasonable numerical 
effort. In typical calculations ts « 1.1 —1.2 fm, while sep¬ 
arations as large as 1.4 fm have been reported only in very 
few cases [121 US] ■ Hence, in order to guarantee a reliable 
determination of Ge and Gm using the plateau method, 
the contributions from excited states in eq. (22) must al¬ 
ready be sufficiently suppressed for t, (ts ~ t) < 0.5 fm. 
Moreover, since the gaps A and A' are proportional to 
in the chiral regime, one expects that this effect will 
become even more pronounced for the more chiral ensem¬ 
bles. 

In Fig.[^ we show effective mass plots for a nucleon 
at rest, computed at two different values of the lattice 
spacing at nearly fixed pion mass. One clearly sees 
that the ground state is isolated only for separations 
larger than 0.5 fm. Since the asymptotic behaviour must 
be reached for both the initial and final-state nucleons, 
which may also carry momentum, source-sink separations 
of the order of 1 — 1.5 fm seem rather small. Therefore, 
one cannot rule out a systematic bias, unless source-sink 
separations significantly larger than I fm are realized. 

In order to minimize or eliminate such a bias in our 
final results, we employ three different methods: 


• Plateau fits: For a fixed value of G the quantities 
t, ts) are fitted to a constant over a small 
interval in t. The default value of the source-sink 
separation is ts « I.lfm. For the high-statistics 
run on the N6-ensemble, we have also considered 
separations as large as ts = 1.4 fm. 


The standard “plateau method” for extracting Ge and 
Gm proceeds by fitting the ratios defined in eq. (13) to 
a constant in the region where they are approximately 
independent of t and G, assuming that their asymptotic 
behaviour has been reached. For the following discussion 
it is useful to define an “effective” electric form factor, 
G^{Q'^,t,ts), by dividing out the kinematical factor in 
eq. ([T4|, i.e. 


GSf = J - -f—Rvoiq,t,ts). (21) 

y rriN + Eq 

A similar relation is used to define G^{Q^,t,ts). As 
t, (tg — t) —> oo, the effective form factors will approach 
their asymptotic values with exponentially small correc¬ 
tions, 

Gfy,{Q\ t, ts) = Ge.m(Q") + G(e-^*) + G(e-^'(‘=-*)), 

( 22 ) 

where A and A' denote the energy gaps between the 
ground and first excited states for the initial- and final- 
state nucleons, respectively. Here we omit the superscript 


• Two-state fits: In this case the leading contribu¬ 
tions from excited states are included explicitly by 
using an ansatz of the form 

Gf^MiQ^,t,ts) = GE,MiQ^) 

+ + Ce.L(Q^) ^ 

with simultaneous fits in t and ts — t performed 
to the data collected for several source-sink sepa¬ 
rations ts- In order to stabilize the fits and reduce 
the number of fit parameters, we fix the gaps to 
A = m,r and A' = 2m,n., assuming that the lowest- 
lying excitations are described by multi-particle 
states, consisting of a nucleon and at least one pion. 
In our chosen kinematics the nucleon at ts is at rest, 
such that the lowest-lying multi-particle state con¬ 
sists of one nucleon and two pions in an S-wave. By 
contrast, the initial state carries momentum and, in 
the absence of irN interactions, therefore consists 
of a moving nucleon and a pion at rest, hence the 
choice A = m,^.. With these assumptions we may 
determine the form factors Ge,m as well as the coef- 
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FIG. 2. Effective masses in the nucleon channel computed at two different lattice spacings for rrijr ~ 330 MeV (left panel) and 
« 275MeV (right panel). 




FIG. 3. Data for G|f computed for several values of ta at the lowest non-zero momentum transfer on the N6 ensemble. The 
yellow band shows the result for Ge determined from the summation method. The solid curves are a representation of the data 
at individual values of ts, as determined from a simultaneous two-state fit to the solid points in the left panel. The grey band 
denotes the corresponding asymptotic value. For the explanation of the dashed curves, see text. 


ficients ^ and as fit parameters for a given 
value of Q^. 

• Summed insertions (“summation method”): Fol¬ 
lowing refs. [23H5H1IM] and our previous work 
EilES] we define the quantities S-e by 

ts—a 

‘S'E,M(Q^^s) := a X! (24) 

t—a 

whose asymptotic behaviour is given by 

-S'e.m(Q^, ts) Ke,m{Q^) + ts G'e,m(Q^) + •. ■, 

(25) 

where Ke,m{Q'^) denote (in general divergent) con¬ 
stants, and the ellipses stand for exponentially sup¬ 
pressed corrections. The precise form of the latter 


depends on the details of the spectrum. If, for in¬ 
stance, A = rriT^ and A' = 2m,n., the leading correc¬ 
tion is of the order exp{—Atg}, while for A = A' 
it is of the generic form 

(^E,M + BE,Mts) exp{—Afg}, (26) 

with coefficients Ae,m and i?E.M- By computing 
SE,Mits) for several sufficiently large values of ts, 
form factors can be determined from the slope of a 
linear fit. Since tg is, by design, larger than either t 
or (ts — t) , excited-state contributions are paramet¬ 
rically reduced compared to the plateau method. 
The summation method has been successfully ap¬ 
plied in our earlier calculation of the nucleon axial 
charge m and also in recent studies of various nu¬ 
cleon matrix elements [2S1 [55] . 
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As a common feature among lattice calculations, we 
note that nucleon electromagnetic form factors are typ¬ 
ically overestimated at a given value of relative to 
the phenomenological representation of the experimen¬ 
tal data m, even when the calculation is performed for 
small pion masses. The three methods which we employ 
to determine G'e,m(Q^) are compared in Fig.|^ Our data 
computed for different source-sink separations ts show a 
systematic downward trend as tg is increased from 0.65 
to 1.1 fm. This reinforces our concern that tg ~ 1.1 fm 
is insufficient to rule out a systematic bias when the 
plateau method is applied. In comparison, the slope de¬ 
termined from the summed ratio in eq. ( |25[ ) yields a result 
for Ge(Q^) which lies sufficiently below tg) for 

all t, tg. 

The result from the two-state fit is even smaller: The 
asymptotic value of G^{Q^,t,ts) for t, {tg — t) —)• oo is 
represented by the grey band, while the dashed curves 
correspond to Ge + c^^ 

where Ge, ^ and ^ are determined from the fit. At 
face value, the sizeable gap between the result from the 
two-state ht and the data for G|f suggest that the latter 
are far from the asymptotic behaviour when tg < 1.1 fm. 
In particular, there is no overlap between the grey band 
and any of the data points from which it is determined. 
The dashed lines in the plot suggest that the two-state 
fit constrains Ge merely from the curvature in t at a 
given tg and from the trend in the source-sink separation 
as the latter is increased. To investigate this further we 
have added two more values of tg to the N6 ensemble, 
corresponding to separations of 1.3 and 1.4 fm, respec¬ 
tively. The additional data for G|f are shown in the 
right panel of Fig.[^ In spite of the large statistical er¬ 
ror, it is clear that G|f approaches the asymptotic value 
extracted from the two-state fit. We conclude that two- 
state fits applied to our data collected for tg < 1.1 fm 
should not simply be discarded, even though the fit cor¬ 
responding to the grey band in the left panel of Fig. id oes 
not appear very convincing. We will thus include such re¬ 
sults in our subsequent analysis, but interpret them with 
the necessary amount of caution. Consequently, our pre¬ 
ferred method for determining form factors remains the 
summation method. 


IV. Q^-DEPENDENCE AND CHIRAL 
BEHAVIOUR 

In this section we discuss the dependence of form fac¬ 
tors on the squared momentum transfer, Q^, and their 
behaviour as the pion mass is tuned towards its physical 
value. Here we focus on the more qualitative features and 
defer a detailed discussion of chiral extrapolations based 
on baryonic ChPT to section|Vj 

For the remainder of this paper, we concentrate on 
results obtained using the conserved (point-split) vector 
current, noting that the local vector current yields fully 
consistent results, provided that it is properly renormal¬ 


ized. Moreover, the ratio of matrix elements computed 
using the local and conserved currents provides an esti¬ 
mate of the renormalization factor Zy, which we find to 
be in agreement with other work m- 


Run 

method 

{rl) [fm^l 

(r-L) [fm^] 

M 

A3 

plat 

0.310(14) 

0.355(29) 

4.22(17) 


sum 

0.335(27) 

0.391(56) 

4.65(34) 


two-state 

0.339(18) 

0.425(55) 

4.48(29) 

A4 

plat 

0.362(23) 

0.324(52) 

3.65(29) 


sum 

0.462(55) 

0.190(63) 

3.15(38) 


two-state 

0.453(35) 

0.216(71) 

3.42(43) 

A5 

plat 

0.413(26) 

0.395(59) 

3.69(31) 


sum 

0.504(57) 

0.333(98) 

3.72(51) 


two-state 

0.543(56) 

0.53(21) 

4.3(1.0) 

B6 

plat 

0.427(22) 

0.442(32) 

3.89(17) 


sum 

0.581(89) 

0.42(10) 

4.11(49) 


two-state 

0.585(56) 

0.69(18) 

4.90(62) 

E5 

plat 

0.304(14) 

0.318(44) 

4.03(30) 


sum 

0.336(22) 

0.373(69) 

4.21(45) 


two-state 

0.385(19) 

0.367(69) 

4.16(48) 

F6 

plat 

0.407(17) 

0.387(27) 

3.67(18) 


sum 

0.451(23) 

0.502(49) 

4.15(25) 


two-state 

0.505(21) 

0.366(68) 

3.32(33) 

F7 

plat 

0.421(25) 

0.385(45) 

3.72(25) 


sum 

0.431(29) 

0.446(62) 

4.17(31) 


two-state 

0.518(35) 

0.54(24) 

3.57(80) 

G8 

plat 

0.463(25) 

0.505(58) 

4.08(30) 


sum 

0.502(76) 

0.45(13) 

4.58(73) 


two-state 

0.739(76) 

0.85(39) 

5.8(1.6) 

N5 

plat 

0.333(12) 

0.314(24) 

3.85(16) 


sum 

0.383(29) 

0.361(60) 

3.93(35) 


two-state 

0.381(12) 

0.399(34) 

4.18(21) 

N6 

plat 

0.391(16) 

0.412(33) 

3.91(21) 


sum 

0.448(19) 

0.336(33) 

3.88(21) 


two-state 

0.571(22) 

0.535(91) 

4.11(43) 

07 

plat 

0.396(19) 

0.399(37) 

3.45(19) 


sum 

0.460(23) 

0.436(42) 

3.66(21) 


two-state 

0.672(38) 

0.99(20) 

4.98(69) 


TABLE II. Electric and magnetic charge radii and magnetic 
moment as determined from the three methods on each of our 
ensembles. 


A full set of results for Ge{Q^) and Gm{Q^) obtained 
from all three methods (i.e. plateau fits, summation 
method and two-state fits) is presented in Tables Vl||XVI 


in appendix]^ In order to describe the dependence on Q" 
we follow the standard procedure of fitting the results for 
Ge,m using a dipole ansatz motivated by vector meson 






























dominance, i.e. 


G.«3»)=(l + g) . 

2-2 (27) 

Gm(<3^) = Gm(0) • 

Using the definition in eq. Q, the charge radii are then 
obtained from 


{rku) = 


12 


^E,M 


(28) 


The ratio Gm(Q^)/Ge(( 5^) is related to the nucleon’s 
magnetic moment /i via 


= 1 + K = 


GuiQk 

Ge(Q2) 


(29) 


Q2=0 


where k denotes the anomalous magnetic moment. The 
reciprocal ratio Ge((5^)/Gm(Q^) is an interesting quan¬ 
tity regarding the discrepancy between experimental de¬ 
terminations based on Rosenbluth separation and the re¬ 
coil polarization technique. 

Examples of dipole fits to the form factor data obtained 
for two different pion masses are shown in Figs.|^ an di 
where they are compared to Kelly’s phenomenological pa¬ 
rameterization |57j of experimental dataj^ In comparing 
experimental and lattice results, one must bear in mind 
that the latter have been obtained at unphysical values 
of the pion mass. 

Clearly, the slope of the electric form factor near 
= 0 varies depending on the method which is used to 
determine Ge from the ratio of correlators. One observes 
that two-state fits produce by far the steepest drop-off, 
while the standard plateau method yields the flattest be¬ 
haviour in Q^. This translates into a corresponding hi¬ 
erarchy for estimates of the electric charge radius, which 
are tabulated in Table |TTj As the pion mass is lowered 
towards its physical value, one also finds that the spread 
in the results for (r^) becomes more pronounced among 
the three methods. This is consistent with the assertion 
that the issue of unsuppressed excited state contributions 
becomes increasingly important near the physical pion 
mass. 

For Gm the systematic trend in the Q^-dependence is 
not so clearly visible as in the case of the electric form 
factor and charge radius, which is partly due to the larger 
statistical errors. Qualitatively, one observes that our lat¬ 
tice data for Gm{Q^) show better overall agreement with 


^ We employ Kelly’s parameterization as a benchmark, since the 
differences between Kelly and the more recent parameterization 
by Arrington and Sick [58], as well as the dispersive analysis by 
Lorenz et al. ESI are too small to be resolved at the level of 
statistical precision provided by our data. 


the representation of the experimental data, regardless of 
the method which the former have been obtained with. 

Dipole fits to the data for Gm(< 3^) extracted from two- 
state fits show a slight - albeit statistically insignificant - 
tendency for larger intercepts at vanishing Q^, resulting 
in somewhat higher estimates for the magnetic moment, 
/r. Moreover, these fits reveal that the electric and mag¬ 
netic radii turn out to be rather similar (see Table[TT| . 
The fact that Gm/Ge shows no statistically significant 
deviation from a constant within the Q^-range we are 
able to investigate is consistent with the experimental 
data extracted using the traditional Rosenbluth separa¬ 
tion technique. 

In order to further assess the effectiveness of the three 
methods employed to extract the form factors, we have 
studied their chiral behaviour at a reference value of the 
four-momentum transfer of = 0.1 GeV^. This is very 

close to the smallest non-zero value of on the G8 en¬ 
semble, i.e. at our smallest pion mass. Using the dipole 
fit parameters, we have obtained GE(Qref) GuiQki) 
for all ensembles in our set. Similarly, we employed the 
phenomenological parameterization to produce the cor¬ 
responding estimates from experiment at and at the 
physical pion mass. The results are shown in Fig.[^ The 
plateau method clearly overestimates Ge, as there is no 
observable tendency for the data at different pion masses 
to approach the experimental result. The summation 
performs slightly better but does not improve the situa¬ 
tion substantially. Only the data based on two-state fits 
show a trend which brings them into agreement with ex¬ 
periment at the physical pion mass. For the magnetic 
form factor the situation is more favourable: the plateau 
method only slightly underestimates Gm relative to ex¬ 
periment, while the chiral trend in the data extracted us¬ 
ing the summation method or two-state fits agrees well 
with Kelly’s parameterization. 

We conclude that the summation method cannot fully 
reconcile lattice data for Ge with its phenomenological 
value. Such an agreement can only be reached if one is 
willing to trust two-state fits. 


V. CHIRAL FITS 

Our task is now to make contact between lattice data 
for form factors obtained for a range of pion masses and 
lattice spacings and the quantities which describe key 
properties of the nucleon, namely the charge radii and 
the magnetic moments. This link is provided by chiral 
effective field theory. The approach which has so far been 
most widely applied to perform chiral extrapolations of 
lattice results for these quantities is based on Heavy 
Baryon Chiral Perturbation Theory (HBChPT) [5D], sup¬ 
plemented by the inclusion of the A-resonance [SIJIS2] • 
Here we employ an alternative formalism, i.e. the man¬ 
ifestly Lorentz-invariant version of baryonic ChPT [29l- 
131], which has also been extended to include the A- 
resonance [53] , as well as vector mesons [311IS3] ■ Our pro- 
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N6 G8 




FIG. 4. Dipole fits of the dependence of Ge, as determined using the plateau method (shown in blue), summed insertions 
(red) and two-state fits (green). The left and right panels correspond to pion masses of 331 and 193 MeV, respectively. The 
black line denotes Kelly’s parameterization of experimental data. 
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G8 



FIG. 5. Dipole fits of the dependence of Gm. 

cedure resembles the strategy pursued in [35] to extract 
charge radii and magnetic moments from experimental 
data of nucleon form factors. In particular, we focus on 
fitting the dependence of the form factors Ge and Gm on 
the pion mass and the squared momentum transfer to 
the expressions of baryonic effective field theory (EFT), 
including vector degrees of freedom. The relevant EFT 
expressions for Ge and Gm have been supplemented by 
terms which describe the dependence on the lattice spac¬ 
ing a. In this way we combine a simultaneous chiral and 
continuum extrapolation with a fit to the Q^-dependence 
of form factors. In order to enable a comparison with the 
standard approach we also perform fits to the pion mass 
dependence of charge radii and the magnetic moment to 


The meaning of the symbols is identical to Fig.[^ 

several variants of HBChPT. 

Specifically we consider the manifestly Lorentz- 
invariant effective Lagrangian describing "kN interactions 
including vector mesons at 0{q^) in the chiral expansion. 
A detailed discussion of this effective Lagrangian can be 
found in |35| . Tablejlllj gives an overview of the various 
interaction terms, as well as the associated low-energy 
constants and hadron masses. From the table one can 
read off which low-energy constants are determined by 
our fitting procedure and which phenomenological infor¬ 
mation is used to fix the values of the remaining ones. 
We note that the interaction terms proportional to c^ 
and dy do not contribute in the iso-vector case consid¬ 
ered here. We have also dropped the contributions from 
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FIG. 6. The Q^-dependence of Gm/Ge- 




ral [GeV^] ral [GeV^J 

FIG. 7. Gomparison of the three methods for extracting the nucleon form factors Ge (left) and Gm (right) at a fixed value of 
. Shown are, from top to bottom, the results from the plateau method, summation method, and two-state fit, as a function 
of m^. 


the tLi-meson entirely, since they were found to have only 
a negligible effect on the results [35]. Furthermore, in 
ref. [65] it was shown that the universal p-meson coupling 
constant g can be fixed via the Kawarabayashi-Suzuki- 
Riadzuddin-Fayyazuddin (KSRF) relation [^|37[, which 
follows by requiring the self-consistency of an effective 
chiral theory involving pion, nucleons and the p-meson. 

The full expressions for the chiral expansions of the 
Dirac and Pauli form factors to 0{q^) are given in ap¬ 
pendix D.2 of ref. [68] but are too lengthy to be displayed 
here. Starting from those formulae, we have formed the 
appropriate linear combinations for the iso-vector form 
factors Ge and Gm- The resulting expressions were 
used to perform a simultaneous fit to both Ge(Q^) and 
Gm(Q^) obtained for a range of pion masses and mo¬ 
mentum transfers, at all three values of the lattice spac¬ 


ing|f] Cutoff effects can be easily incorporated into this 
framework by adding terms proportional to the lattice 
spacing a to the form factors, i.e. 

Ge(Q') = Gi^^(Q2) + aQ2^E, 

Gm(Q") = G^'^T(Q2) + apM, (30) 

where G|^J denote the continuum EFT expressions for 
the form factors, while the coefficients /3e,m are taken as 
fit parameters. This ansatz takes account of the fact that 
the matrix element corresponding to the electric form 
factor is G(a)-improved at vanishing momentum transfer. 


® To evaluate the loop integrals appearing in the expressions, we 
make use of LoopTools 16911701 . 
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Interaction 

Low-energy 

parameter 

Value 


F 

= 92.2 MeV 


Ml 

Lattice input 


rhp 

^exp _ yyg MeV or lattice input 


9 

g = mp/'/2F-n = 5.93 from KSRF relation 


dx 

Fit parameter 

^7TpN 

9 

from KSRF relation 


Gp 

Fit parameter 

r(^) 

^7tN 

rfiN 

— 93 g MeV or lattice input 

5A 

9T = 1-27 

r(2) 

^7tN 

C6 

Fit parameter 


C7 

C 7 does not contribute in the iso-vector case 

r(^) 

^7tN 

da 

Fit parameter 


dj 

d-j does not contribute in the iso-vector case 

■ TIT) 

^ttA 

TTlA 

= 1210 MeV 

<1 

w b: 

QttN 

1.125, from fit to A — >■ nN decay width 


TABLE 111. Interaction terms of the chiral effective theory used for fitting lattice data for nucleon electromagnetic form factors. 
A circle on top a symbol indicates that the corresponding low-energy parameter is defined in the chiral limit. The third column 
describes how their values are determined. Further details are described in the text. 


Estimates for the charge radii (4), and the 

anomalous magnetic moment k are obtained by insert¬ 
ing the fitted values of the low-energy parameters of 
dx^Gp^d^ and cq = cq — Gpj^g into the corresponding 
EFT expressions, i.e. 


/ 2 \ __ 6 

k = Gm(0)-1. 


dGE.M(g') 

dQ^ 


Q2=0 


(31) 

(32) 


The relations between these quantities and the Dirac ra¬ 
dius (rj), as well as the combination ^(r^) reads 


(-?> = (4 

«(^2) = 


^ 4to^’ 

(1 -b k) (r^) - (r|) -b 


6k 

4m^ 


(33) 


We refer to a fit applied to form factor data obtained 
from the summation method, over the entire range of 
with a pion mass cut of = 330 MeV, and the 

masses of the p-meson and nucleon fixed to their exper¬ 
imental values as our standard procedure. Standard fits 
were performed with and without terms parameterizing 
lattice artefacts (see eq. (HI). For the fit including lattice 
artefacts one finds Xred = ^-21, for 66 degrees of freedom. 
Results for the quantities (r|), (r^), k, {rf) and the com¬ 
bination /c(r'l) are listed in Table [l\^ One observes that 
fits with and without 0{a) terms produce compatible re¬ 
sults: differences at the level of at most 1.5 standard 
deviations are seen in (r^) and k. This indicates that 
the influence of lattice artefacts on the results is small. 

A first qualitative assessment can be made by plot¬ 
ting the Q^-dependence of Ge and Gm at the physical 


point compared to various phenomenological parameter- 
izations of experimental data, as shown in Fig.[^ While 
the dependence of Ge on the squared momentum trans¬ 
fer is somewhat flatter compared to experiment, the be¬ 
haviour of Gm is reproduced very well. 

In order to estimate the systematic error, we have con¬ 
sidered a number of variations in the fitting procedure, 
which are compiled and labelled in Table[V| These in¬ 
clude different pion mass cuts, restrictions of the fitted 
range in and the use of the masses of the p and nucleon 
determined by the lattice calculation at the respective 
value of the pion mass. These variations are indicative of 
higher-order terms in the chiral expansion and probe the 
overall consistency of our particular EFT approach. We 
have also estimated the residual systematic uncertainty 
due to excited states, by repeating the entire procedure 
using the form factor data obtained from two-state fits. 
Variations of the fitting proced ure (corresponding to the 
entries in lines 3 — 7 in Table [lv| were always applied 
neglecting terms parameterizing lattice artefacts (i.e. for 
/3 e = /3m = 0), as this produced more stable fits, in par¬ 
ticular when imposing more aggressive cuts in the pion 
mass or range. However, while the systematic error 
budget is estimated from fits excluding lattice artefacts, 
we prefer to quote our main results using fits in which 
G(a)-terms have been accounted for. 

We thus obtain as our final results: 

(r|) = 0.722 ± 0.034 (stat);°:°30 (xfit)t°:°®g (exc) fm", 
(4) = 0.720 ± 0.053 (stat)tn25 (exc) fm^, 

K = 3.33 ± 0.35 (stat);^;4^ (xfit)io.?6 (exc), (34) 

(rl) = 0.501 ± 0.041 (stat)t;);(!^^ (xfit)to:ooo (exc) W, 

K (r|) = 2.61 ± 0.09 (stat))/^);?^ (xfit)io(o5 (exc) fm^. 
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FIG. 8. The dependence of Ge and Gm at the physical pion mass and in the continuum limit, as determined from a 
simultaneous fit to lattice data (“direct” fit) including lattice artefacts. The bands denote the statistical error. The solid, 
dashed and dashed-dotted curves are the phenomenological representations of experimental data of refs. [57],[Ml and | 59| . 
respectively. 


(^E)[fm1 





Fit 

0.722(34) 0.720(53) 3.33(35) 

0.501(41) 

2.61(9) 

standard, with G(a)-terms 

0.748(12) 

0.636(8) 

3.93(11) 

0.487(14) 

2.65(9) 

standard, without G(a)-terms 

-0.013 

-0.004 

-0.15 

-0.003 

-0.11 

variant 1 

+0.030 

-0.025 

+0.28 

+0.011 

+0.04 

variant 2 

-0.001 

-0.009 

-0.14 

+0.009 

-0.14 

variant 3 

+0.030 

+0.035 

-0.40 

+0.057 

-0.16 

variant 4 

+0.083 

+0.029 

-0.16 

+0.093 

-0.05 

variant 5 


TABLE IV. Results for charge radii and magnetic moments from direct fits to the form factors. Lines 3 — 7 contain the 
differences between variants of the fitting procedure as labelled in TablefVj and the results from the standard fit performed 
without C)(a)-terms. For instance, variant 5 yields {r^) = 0.831. 


Here, the systematic uncertainties estimated from fit 
variants 1 — 4 have been combined into an overall chi¬ 
ral fitting error, while the difference between employing 
the summation method and two-state fits is quoted as 
a separate, residual systematic uncertainty arising from 
excited states. 

We did not consider hts to baryonic EFT including 
A degrees of freedom when assessing our systematic er¬ 
rors, as such fits produced unacceptably large values of 
Xied ’'^hen the low-energy parameter g-^NA was fixed to 
the phenomenological value of 1.125. On the other hand, 
treating gT^NA as a ht parameter resulted in an unphysi- 
cally small value. 

In Figs.j^andlO we compare the estimates for (r|), 
(r^) and k at the physical point (shown as yellow points) 
with experiment. While and k agree quite well with 
the experimental results within statistical errors, we find 
that direct fits to the form factors underestimate the elec¬ 
tric radius. However, given the large systematic uncer¬ 
tainty, we note that our estimate for (r|.) is not incom¬ 
patible with either the CODATA result nni or the value 


determined from muonic hydrogen [HlISj. 

In order to enable a comparison with previous lattice 
determinations [nH23], we have performed chiral extrap¬ 
olations of our data to the expressions of HBChPT in¬ 
cluding the A-resonance at 0{e^) in the small scale ex¬ 
pansion m- In particular, we fitted the pion mass depen¬ 
dence of our data for the Dirac radius (rf), the anoma¬ 
lous magnetic moment n and the combination to 

the expressions of refs. [SI isi mi, which are summa¬ 
rized, e.g. in appendix A of[23j. In these fits the low- 
energy parameters gA, and g-^NA have been fixed to 
the same values as in Table|Tnj An additional parameter, 
the photon-nucleon-A coupling, was fixed to the value 

-2.26 [S]. 

The results of such a HBChPT fit, with a pion mass 
cut of 330 MeV, without terms parameterizing cutoff ef¬ 
fects and with the masses of the nucleon and A fixed 
to mN = 938 MeV and ttia = 1210 MeV, respectively, 
are shown in Fig.|^ While the value for (r^) agrees 
within statistical errors with the result determined from 
directly fitting the form factors, there is a deviation by 
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Standard fit 

Variants 

Label 

Impose pion mass cut of < 330 MeV 

mT < 300 MeV 

1 


No mass cut 

2 

Fit entire available range in Q'‘ 

Impose cut of < 0.5 GeV'' 

3 

Use experimental values for mti,mp 

Use lattice input for m^,mp 

4 

Fit data obtained using the summation method 

Fit data extracted from two-state fits 

5 


TABLE V. The standard procedure for fitting nucleon form factors to the expressions from baryonic ChPT and the variants 
applied in order to estimate the systematic error. 




FIG. 9. Pion mass dependence of electric and magnetic radii extracted by fitting the Q^-dependence of form factors to a 
dipole form. Chiral fits to HBChPT for m,r < 330 MeV and their statistical uncertainty are represented by the bands. The 
yellow points denote the results obtained by directly htting the form factors to the EFT expressions. The associated systematic 
uncertainties are shown by the dashed error bars. 



ml [GeV"] 

FIG. 10. Pion mass dependence of the anomalous magnetic moment. For an explanation of symbols, see Fig.|^ 


more than two standard deviations in the case of (rg), 
which, however, becomes insignificant when systematic 
errors are taken into account. Interestingly, the result 
for the electric radius obtained from the HBChPT fit is 
statistically compatible with both the CODATA estimate 
and the value determined from muonic hydrogen. How¬ 
ever, a number of comments are in order: hrstly, we note 


that HBChPT fits including terms parameterizing lattice 
artefacts mostly fail. This may be due to the lack of a 
clear trend in the data for the charge radii as the lat¬ 
tice spacing is varied. Secondly, hts based on HBChPT 
depend much more strongly on whether the input data 
originate from applying the summation method or two- 
state hts. 
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We note in passing that we have also applied baryonic 
EFT based on the Lagrangian described in Table [TlT| to 
perform chiral extrapolations of charge radii and n as 
an alternative to HBChPT. The results at the physical 
point are consistent with direct fits to the form factors 
within statistical errors, except for the anomalous mag¬ 
netic moment. We conclude that the chiral behaviour 
of the quantities computed here is not fully understood 
in terms of baryonic chiral effective theory. In order to 
clarify the situation it is mandatory to obtain more sta¬ 
tistically precise data at the physical pion mass. 

Given that the results for the charge radii and k ex¬ 
hibit a large spread depending on the details of the chiral 
fitting procedure, the agreement of the HBChPT result 
for (r|,) with the experimental values must be considered 
an accident. Due to the better overall stability we prefer 
to quote our final estimates from fits applied directly to 
the form factors Ge and Gm (see eq. (34)). 


VI. CONCLUSIONS AND OUTLOOK 


We have performed a comprehensive study of the iso¬ 
vector electromagnetic form factors of the nucleon in two- 
flavour QCD with a strong emphasis on controlling the 
various sources of systematic error. Our findings culmi¬ 
nate in the estimates shown in eq. (34) for the charge radii 


and magnetic moment including a full error budget. 

While the summation method provides a handle to 
explore excited-state contributions independently of the 
standard plateau method, the issue of a systematic bias 
could not be fully resolved. Although we prefer the sum¬ 
mation method, since it does not rely on specific values 
of the energy gaps, two-state fits produce values that 
appear to reproduce the phenomenological situation bet¬ 
ter. However, for lack of data at source-sink separations 
of more than 1.5 fm, we cannot currently resolve the is¬ 
sue completely. Still, our data support the notion that 
agreement with experiment can be obtained by mapping 
out the pion mass dependence close to the physical point 
in conjunction with addressing the issue of excited-state 
contaminations, which has a greater impact as the pion 
mass is reduced. 

For the first time, we publish a complete error budget 
for baryonic charge radii and magnetic moments. Also 
for the first time, we have applied the method of [3^, i.e. 
applying the full framework of covariant baryonic chiral 
effective theory to the nucleon electromagnetic form fac¬ 
tors, in the context of lattice QCD. 

While the wider picture of the convergence properties 
of the various forms of baryonic chiral effective theory 
cannot be fully addressed with our data, we have a clear 
preference for applying the fully covariant chiral effective 
theory to the form factors themselves. In particular, this 
avoids the systematic uncertainties inherent in the use of 
the somewhat simplistic dipole description of the form 
factors. 

In order to further resolve the issue of excited-state 


effects, significant improvements in statistical precision 
at larger source-sink separations will be required. One 
proposed way to achieve this is the use of methods such 
as All-Mode Averaging (AMA) [75] in order to reduce 
the associated computational costs. First studies have 
yielded encouraging results, and we intend to pursue this 
avenue further. Another direction for improvement will 
be to make use of the variational method [5il I75H75] for 
increased control and suppression of excited-state contri¬ 
butions. 

Going beyond the iso-vector form factors to the pro¬ 
ton and neutron form factors will require the inclusion of 
quark-disconnected diagrams, which makes it necessary 
to use all-to-all propagators. By employing stochastic 
estimators along with a generalized hopping parameter 
expansion (gHPE), the scalar charge radius of the pion, 
which receives a significant contribution from quark- 
disconnected diagrams, has recently been measured on 
the lattice with an overall accuracy rivalling that of phe¬ 
nomenological determinations [TiH - ITH] . Related methods 
have been employed by other groups to study the nucleon 
structure quantities with disconnected parts iiniiMsil, 
and we intend to further address the proton radius puzzle 
from the lattice side by using these and similar methods 
to separately study proton and neutron form factors in 
the future. 
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Appendix A: On the use of non-covariant 
interpolating operators in nucleon form factor 
calculations 


Typically, the interpolating nucleon operators used in 
lattice form-factor calculations are not Lorentz covariant, 
due to the spatial smoothing procedures applied on the 
elementary fields. If we denote the overlap of the inter¬ 
polating operator 'I'q)^) onto the nucleon as followt]^ 

(0|vI/„(0,a;)|A,p,s) = C/^(p)e*P'", (Al) 

the question then poses itself, whether the fact that 
U^{p) does not obey the Dirac equation affects the cal¬ 
culation in any way. 

We denote by u^(p) the usual plane-wave solutions to 
the Dirac equation^ normalized according to 


3 Ep xq 

2E„ 


Y^KipMip), (A5) 




A 


g EpXQ 

2E„ 


(A6) 


C^2.isiP^^o) = j d^xe ('0a(a;)4'/3(O)) 

o EpXQ _ ^ 

Y^KipWliP)- (A7) 


= z. 


We now define 


Mip) = '^U%p)u%p). 


(A8) 


The other spin sums can also be expressed through Ai. 
First, we have 

^u%p)iE{p) = 7oA^(p)'''7o- (A9) 


u'^{p)u^{p) = 2mAr(5’’® 
and obeying the spin sum rule 

u%p)u^{p) = Epjo - ip ■‘y + niN- 

s^l,2 


(A2) Since the overlap of the local nucleon operator is a gen¬ 
uine Dirac spinor, the matrix M satisfies 


M{p){joEp - ip ■'f - rriN) = 0- (AIO) 

(A3) 

We observe from the definition of A4 that 


Let ipaix) be a nucleon interpolating field which does 
transform as a covariant Dirac spinor. Then its overlap 
onto the nucleon state has the form 

{0\tpa{0,x)\N,p,s) = Ziul^{p)EP'^, (A4) 


M{p)u%p) = 2mNU%p), 


and hence also 


U%p) 


I 

2m iq 


u^phoMipY-fo- 


(All) 


(A12) 


where Zi is independent of p and can be chosen real and 
positive by an appropriate choice of the phase of the field 
il)a{x). To answer the question formulated above, we 
analyze the covariance properties of correlation functions 
of the operators ^'^(a;) and ipaix). We will focus on the 
asymptotic behaviour of the correlation functions at large 
Euclidean time separations, where they are saturated by 
the nucleon ground state, and we indicate this by a * in 
the equations below. 

We consider the following two-point functions and 
their spectral representation for xq > 0 , 

C2,LiP^xo) = J d^x (4'„(a;)^/j(0)) 


We can thus write the spin sum appearing in the two- 
point function C 2 ,ss as 


U\p)U\p) = -^M{p){Ep-iQ - fp • 7 -h toat) 

X 7oA4(p)'f7o (A13) 

= 7^AI(p)7oAI(p)'''7o. (A14) 

2m N 


The second equality uses Eq. (AIO). Thus 70 times the 
spin sum of the C/’s is a hermitian matrix. 

Let J{x) be a local operator with the following matrix 
elements between one-nucleon state^ 


{N,p\s'\J{0)\N,p,s) = u"' ip')Jiq)u^p), (A15) 


^ The states are normalized according to (N,p', s'\N,p, s) = 
2Ep5,,, (27r)35(3)(p'-p). 

® In this appendix we use the conventions of Peskin and Schroeder 
in Minkowski space (with Dirac matrices ) and introduce a 
corresponding set of Euclidean Dirac matrices, 70 = 7 m 
7 fc = The are all hermitian; 70 , 72 and 75 are sym¬ 

metric, 7 i and 73 are antisymmetric. 


where q = p' — p. The three-point function relevant to 
form factor calculations with vanishing momentum at the 


For instance, J{q) = = 

IK ’ 7 m]) vector current i/>(fc) 7 j^ ip(x). 
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sink and its spectral representation read (for xq > yo > 

0 ) 


C 3 ^j{Q,yo,xo) = J (fy J <fx J{y) ^'/3(0)) 

■(5]C7«'(0)<(0) 


^-EqVo g-m(£Co-i/o) 


2Eq 2m ff 


X {JiQ))^s{'l2'^si-Q)Upi-q)y (A16) 


The three-point function, projected with a generic pro¬ 
jector r, 

p-EqVo p-m(xo-yo) 

= (A17) 

X Tt{rA<(0)J-(9)7„A<(-9)t7„} 

can thus also be expressed in terms of the matrix At (p). 
The latter can be expanded in the 16 linearly independent 
spinor-space matrices, and symmetries can be used to 
restrict the terms that contribute. 


1. Symmetry constraints on nucleon two-point 
functions 


Let 

C 2 , 0 x(p> 2 ;o) = J (fix {<j){x)x{0)) (A 18 ) 

be a generic nucleon two-point function with interpolat¬ 
ing operators (p and y. We assume that the latter are 
good spinors with respect to spatial rotations and with 
respect to all discrete symmetries, but no assumption is 
made on their behaviour under boosts. The implications 
of the various symmetries are: 

• rotation invariance: 

C2,4,x(Rp, a::o) = Ai {R)C2,^xiP^ ^o)Ai (A)"\ (A19) 

where —>■ Ki{R)(j){R~^x) is the transforma¬ 

tion law of the quark fields; 

• parity: 

C'2.0x(“P>^o) = 7 oC'2,0x(p> 2 : 0)70 ; (A20) 

• Euclidean time reversal: 

C2>x(p> - 2 : 0 ) = 7075 C'2,0x(p> 2 : 0 ) 7570 ; (A21) 

• charge conjugation: assuming that the nucleon in¬ 
terpolating fields transform like the quark fields 
{qa{x) (g(a:)7o72)a, qa{x) (7o72<?(a;))afl 

C'2,0x(-P>-2:o)^ = 7270 C'2.x<^.(P: 2 : 0 ) 7072 - (A22) 


^ This condition fixes the phase of the interpolating fields. An 
interpolating field that satisfies this condition is Xa{x) = 
e“*’'=(iig(707275)^7d!^)7‘a, but e®'<°Xa(a;) does not for (,5 7 ^ 0 , 7 r. 


Combining the three discrete symmetries, we obtain, for 
later use, 

75C'2,,^x(p>2:o)^75 = 727oC2.x0(P: 2 : 0 ) 7072 - (A23) 

The most general form of the nucleon two-point func¬ 
tions allowed by rotation symmetry and parity is 

C 2 , 4 ,xiP,xo) = Rs{p^,xo) -h Jb(p^,a:o)7o (A24) 

-iRv{p'^,xo)p ■ 7 - 2 : 0)70 P ■ 7- 

Time-reversal invariance implies that J-q and J-t are odd 
functions of xq, while Rs and J-y are even functions of cco- 
Charge conjugation relates the functions R parametriz¬ 
ing the correlator C'2,0x those parametrizing the cor¬ 
relator C2,x4>- 


2. General parametrization of Af(p) 


In view of (A24), we parametrize M.{p) as 

M{p) = mNfs{p^) + /o(p^)Ap 7 o - ifv{p^)p ■ 7 

- 2 /t(p^)p-( 707 )- (A25) 


The condition (AlO) from the spectral representation im¬ 
plies the constraints 


fv{p^) = ^(a^/o(p") - rnlfsip^)), (A26) 
/t(p^) = [fs{p^) - Mp^))- (A27) 


In particular, fo{p^) — fs[p^) = O(p^); the case of a 
Lorentz-covariant spinor source corresponds to fo = fs = 
fv independent of p and fr = ^ identica lly. N ow, com¬ 
bining the discrete-symmetry property (A23) and the 
spectral representations (A6 A7), one derives the prop¬ 
erty that /o, /s, fv and /t must all be real-valued func¬ 
tions. 

It is convenient to decompose A4 as follows. 


M{p) = M+{p) + M-{p), (A28) 

M±{p) = \{l±xo)M{p). (A29) 

Using the constraints, one finds the following general 
form of M±{p), 

M + {p) = 5(1 + 70 ) Z+{p){Ep + mN -ip- 7 ), (A30) 
M-(p) = 5(1 - 7 o) Z^{p){Ep - mN + ip-l) (A31) 
with the relations 


^ mNfs{p^) + Epfo[p^) 
Ep + mjx 

^ ^ rriNfsip^) - Epfojp^) 

Ep ruTv 


(A32) 

(A33) 


The bottom line is that the matrix Ai (p) is parametrized 
by two (spatially scalar) functions Z±{p), which are 
linear functions of the nucleon interpolating operator 
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Its phase can be chosen such that Z±{p) are 
real for all p. Furthermore it can be chosen such that 
Z+(p = 0) is positive. The generic case that we will 
consider is then that it remains positive for all mo- 
mentaj^ As a side remark, we note that in the rest frame 
only parity-odd states (not considered here) contribute 
to Tr{(l - 7 o)C2,sz(p = 0,a;o)}. 


3. Final form of the non-covariant two-point and 
three-point functions 


We now specialize to the projector 

r = 5(1 + 7 o)(i + *7573)- 


(A34) 


Starting from Eq. (A5) and _|A13 ), the projected two- 
point function can be writterr] 


Tr{rC2.ss(p,xo)} = 


o EpXQ 


1 


2Ep 2m N 


Tr{(l -f 17573 ) 


X Af + (p)7oAl+(p)^7o} (A35) 
= \Z+{p)f{l + mM/Ep)e-^-^°. 


The second equality uses the parametrization Eq. (A30). 


The term 17573 in Fp does not co ntribu te to this expres¬ 
sion. The trace appearing in Eq. ( |Al7t takes the form 

Tr{FpAf(0)J(g)7oM(-q)^7o} = (A36) 

Tr |(1 -k z 7573 )Af+( 0 )J^(q) 7 oAf+(-q)^|. 

Thus the expression for the three-point function be¬ 
comes 


C3{q,yo,xo) = 


g-EqVo g-mN{xo-yo) 


2Eq 2mjq 


Z+{Q)Z+[-qy 


rriM Tr |(1 -k 7 o)(l -k * 7573 ) J"(<7)(wjv -k Eg -k *q • 7 

(A37) 

As we have seen, the phase of the nucleon interpolating 
operator d* (x) can be chosen such that (p) > 0 for all 
p. If the phase had not been chosen in this way, the 
common phase would nonetheless cancel in the product 
Z+{Q)Z+{-qY. 

For a covariant source, the result would be identical to 


(A37), except that Zj^{p) would be independent of p. In 


the standard expression of the ratio (131, the three-point 


function is however divided by the appropriate combina¬ 
tion of two-point functions to cancel the overlap factor 
Z+{p) for each value of p. The correctness of our cal¬ 
culation is thus not affected by the use of non-covariant 
interpolating operators. 


Appendix B: Tables 


In Tables VI - XVI we give all of our results for the 
iso-vector vector form factors Ge and Gm of the nucleon 
at all values of measured on each ensemble. Listed 
in each case are the values obtained using the plateau 
method with a source-sink separation of tg = 1.1 fm, 
the summation method, and an explicit two-state fit (cf. 
the main text for details). The statistical errors on each 
data point are quoted in parentheses following the central 
value. 


® A priori, it could happen that it becomes negative for some mo¬ 
menta, but then by continuity it would have to vanish some¬ 
where, and there is no symmetry reason for this to happen. 

® In the following equations we allow for the case that the Z+(p) 
have a common phase for all p. 
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A3 

Ge 

Gm 

[GeV^] 

Plat (1.1 fm) 

Summation 

Two-state 

Plat (1.1 fm) 

Summation 

Two-state 

0.230 

0.752 (0.011) 

0.725 (0.018) 

0.734 (0.012) 

3.040 (0.100) 

3.320 (0.175) 

3.067 (0.128) 

0.443 

0.601 (0.013) 

0.583 (0.025) 

0.592 (0.017) 

2.359 (0.078) 

2.503 (0.147) 

2.272 (0.111) 

0.639 

0.495 (0.018) 

0.480 (0.037) 

0.484 (0.022) 

1.921 (0.084) 

1.869 (0.176) 

1.783 (0.125) 

0.823 

0.373 (0.023) 

0.322 (0.056) 

0.369 (0.036) 

1.579 (0.098) 

1.716 (0.216) 

1.536 (0.158) 

0.995 

0.331 (0.023) 

0.299 (0.055) 

0.320 (0.030) 

1.349 (0.091) 

1.416 (0.205) 

1.253 (0.143) 

1.156 

0.283 (0.031) 

0.228 (0.073) 

0.213 (0.037) 

1.194 (0.121) 

0.881 (0.295) 

0.910 (0.162) 


TABLE VI. A3 ensemble (a = 0.079 fm, = 473 MeV): Vector form factors at all values for all extraction methods. 


A4 

Ge 

Gm 

Q" [GeV^] 

Plat (1.1 fm) 

Summation 

Twc 

i-state 

Plat (1.1 fm) 

Summation 

Two 

i-state 

0.229 

0.720 

(0.024) 

0.685 

(0.037) 

0.668 

(0.021) 

2.806 

(0.147) 

2.603 

(0.232) 

2.726 

(0.239) 

0.437 

0.547 

(0.022) 

0.466 

(0.043) 

0.489 

(0.031) 

2.121 

(0.109) 

2.360 

(0.210) 

2.518 

(0.210) 

0.628 

0.459 

(0.024) 

0.380 

(0.048) 

0.385 

(0.045) 

1.936 

(0.130) 

2.166 

(0.199) 

2.146 

(0.246) 

0.805 

0.401 

(0.039) 

0.397 

(0.087) 

0.339 

(0.076) 

1.591 

(0.239) 

1.499 

(0.339) 

1.408 

(0.375) 

0.970 

0.338 

(0.025) 

0.338 

(0.080) 

0.291 

(0.052) 

1.317 

(0.128) 

1.449 

(0.281) 

1.594 

(0.262) 

1.123 

0.299 

(0.035) 

0.242 

(0.099) 

0.205 

(0.077) 

1.046 

(0.156) 

1.192 

(0.352) 

1.510 

(0.409) 


TABLE VII. A4 ensemble (a — 0.079fm, rriw = 364MeV): Vector form factors at all values for all extraction methods. 


A5 

Ge 

Gm 

[GeV^] 

Plat (1.1 fm) 

Summation 

Two 

I-state 

Plat (1.1 fm) 

Summation 

Two 

-state 

0.228 

0.692 

(0.018) 

0.651 

(0.036) 

0.629 

(0.039) 

2.572 

(0.159) 

2.619 

(0.270) 

2.672 

(0.349) 

0.434 

0.521 

(0.019) 

0.461 

(0.038) 

0.433 

(0.040) 

1.967 

(0.121) 

2.082 

(0.203) 

2.021 

(0.272) 

0.623 

0.424 

(0.029) 

0.305 

(0.058) 

0.328 

(0.054) 

1.720 

(0.134) 

2.115 

(0.277) 

1.408 

(0.328) 

0.797 

0.343 

(0.043) 

0.269 

(0.096) 

0.171 

(0.087) 

1.247 

(0.158) 

0.673 

(0.426) 

1.430 

(0.442) 

0.959 

0.255 

(0.033) 

0.176 

(0.084) 

0.002 

(0.075) 

1.023 

(0.135) 

0.989 

(0.350) 

0.867 

(0.381) 

1.110 

0.065 

(0.093) 

0.189 

(0.114) 

0.067 

(0.128) 

0.760 

(0.340) 

0.750 

(0.523) 

-0.338 

(0.522) 


TABLE VIII. A5 ensemble (a = 0.079 fm, 


316MeV): Vector form factors at all values for all extraction methods. 


B6 

Ge 

Gm 

[GeV^] 

Plat (1.1 fm) 

Summation 

Twc 

I-state 

Plat (1.1 fm) 

Summation 

Two 

I-state 

0.104 

0.826 

(0.008) 

0.785 

(0.023) 

0.760 

(0.018) 

3.264 

(0.121) 

3.692 (0.370) 

3.794 

(0.284) 

0.203 

0.709 

(0.012) 

0.664 

(0.036) 

0.669 

(0.025) 

2.819 

(0.121) 

3.172 (0.258) 

2.571 

(0.260) 

0.297 

0.622 

(0.017) 

0.572 

(0.048) 

0.541 

(0.038) 

2.312 

(0.121) 

2.417 (0.245) 

2.283 

(0.299) 

0.387 

0.565 

(0.018) 

0.569 

(0.069) 

0.474 

(0.042) 

2.154 

(0.132) 

2.435 (0.241) 

2.115 

(0.300) 

0.474 

0.507 

(0.018) 

0.460 

(0.053) 

0.428 

(0.044) 

2.034 

(0.105) 

2.149 (0.216) 

1.970 

(0.247) 

0.557 

0.445 

(0.020) 

0.421 

(0.059) 

0.374 

(0.056) 

1.656 

(0.069) 

1.680 (0.208) 

1.399 

(0.233) 


TABLE IX. B6 ensemble (a = 0.079fm, — 268 MeV): Vector form factors at all values for all extraction methods. 
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E5 

Ge 

Gm 

[GeV^] 

Plat (1.1 fm) 

Summation 

Twc 

i-state 

Plat (1.1 fm) 

Summation 

Two 

i-state 

0.356 

0.663 

(0.013) 

0.635 

(0.018) 

0.601 

(0.014) 

2.611 

(0.105) 

2.544 (0.150) 

2.540 

(0.139) 

0.675 

0.477 

(0.016) 

0.443 

(0.025) 

0.411 

(0.018) 

1.909 

(0.090) 

1.773 (0.130) 

1.782 

(0.111) 

0.966 

0.379 

(0.026) 

0.357 

(0.039) 

0.339 

(0.029) 

1.429 

(0.104) 

1.346 (0.161) 

1.350 

(0.151) 

1.233 

0.209 

(0.077) 

0.197 

(0.065) 

0.204 

(0.049) 

1.342 

(0.358) 

1.015 (0.242) 

1.016 

(0.210) 

1.480 

0.220 

(0.046) 

0.191 

(0.060) 

0.190 

(0.038) 

0.996 

(0.178) 

0.887 (0.242) 

0.929 

(0.172) 

1.709 

0.204 

(0.086) 

0.164 

(0.082) 

0.185 

(0.054) 

0.454 

(0.301) 

0.644 (0.301) 

0.876 

(0.241) 


TABLE X. E5 ensemble (a = 0.063 fm, m-^ 


457 MeV): Vector form factors at all values for all extraction methods. 


F6 

Ge 

Gm 

Q" [GeV^] 

Plat (1.1 fm) 

Summation 

Twc 

I-state 

Plat (1.1 fm) 

Summation 

Two 

I-state 

0.162 

0.766 

(0.009) 

0.746 

(0.011) 

0.725 

(0.009) 

2.834 

(0.125) 

2.963 

(0.143) 

2.586 

(0.216) 

0.314 

0.627 

(0.012) 

0.604 

(0.015) 

0.555 

(0.014) 

2.410 

(0.102) 

2.350 

(0.116) 

2.128 

(0.166) 

0.457 

0.541 

(0.016) 

0.516 

(0.019) 

0.455 

(0.020) 

2.028 

(0.103) 

1.809 

(0.128) 

1.962 

(0.173) 

0.593 

0.423 

(0.017) 

0.394 

(0.023) 

0.332 

(0.027) 

1.586 

(0.084) 

1.505 

(0.116) 

1.474 

(0.199) 

0.723 

0.383 

(0.017) 

0.343 

(0.024) 

0.307 

(0.025) 

1.435 

(0.072) 

1.245 

(0.111) 

1.306 

(0.165) 

0.846 

0.342 

(0.021) 

0.291 

(0.028) 

0.253 

(0.030) 

1.278 

(0.076) 

1.007 

(0.133) 

1.205 

(0.171) 


TABLE XL F6 ensemble (a = 0.063fm, = 324MeV): Vector form factors at all values for all extraction methods. 


F7 

Ge 

Gm 

[GeV^] 

Plat (1.1 fm) 

Summation 

Twc 

I-state 

Plat (1.1 fm) 

Summation 

Two 

I-state 

0.162 

0.759 

(0.013) 

0.751 

(0.016) 

0.723 

(0.017) 

2.894 

(0.167) 

3.106 (0.192) 

2.637 

(0.329) 

0.313 

0.623 

(0.020) 

0.612 

(0.020) 

0.566 

(0.025) 

2.383 

(0.136) 

2.462 (0.158) 

1.893 

(0.258) 

0.455 

0.512 

(0.026) 

0.507 

(0.027) 

0.421 

(0.033) 

1.944 

(0.111) 

2.065 (0.158) 

1.613 

(0.294) 

0.589 

0.442 

(0.026) 

0.410 

(0.033) 

0.270 

(0.048) 

1.706 

(0.130) 

1.649 (0.175) 

1.284 

(0.307) 

0.717 

0.417 

(0.037) 

0.372 

(0.031) 

0.289 

(0.048) 

1.461 

(0.107) 

1.465 (0.141) 

1.336 

(0.302) 

0.838 

0.336 

(0.036) 

0.308 

(0.038) 

0.200 

(0.051) 

1.295 

(0.126) 

1.329 (0.165) 

0.911 

(0.342) 


TABLE XII. F7 ensemble (a = 0.063fm, m.^ = 277MeV): Vector form factors at all values for all extraction methods. 


G8 

Ge 

Gm 

[GeV^] 

Plat (1.1 fm) 

Summation 

Two-state 

Plat (1.1 fm) 

Summation 

Two-state 

0.092 

0.870 (0.013) 

0.928 (0.041) 

0.771 (0.030) 

3.367 (0.243) 

4.272 (0.618) 

4.555 (0.815) 

0.181 

0.731 (0.014) 

0.737 (0.046) 

0.563 (0.040) 

2.822 (0.160) 

3.293 (0.398) 

3.081 (0.546) 

0.266 

0.629 (0.016) 

0.591 (0.055) 

0.396 (0.060) 

2.464 (0.155) 

2.972 (0.383) 

2.879 (0.600) 

0.347 

0.558 (0.022) 

0.510 (0.071) 

0.292 (0.074) 

2.173 (0.131) 

3.088 (0.410) 

2.168 (0.568) 

0.426 

0.499 (0.018) 

0.473 (0.058) 

0.149 (0.065) 

1.902 (0.119) 

2.311 (0.315) 

1.862 (0.483) 

0.503 

0.458 (0.020) 

0.419 (0.068) 

0.048 (0.082) 

1.688 (0.130) 

1.899 (0.370) 

1.767 (0.492) 


TABLE XIII. G8 ensemble (a = 0.063fm, m.^- = 193MeV): Vector form factors at all values for all extraction methods. 







20 


N5 

Ge 

Gm 

[GeV^] 

Plat (1.1 fm) 

Summation 

Twc 

i-state 

Plat (1.1 fm) 

Summation 

Two 

i-state 

0.257 

0.708 

(0.008) 

0.678 

(0.018) 

0.687 

(0.008) 

2.858 

(0.086) 

2.703 (0.164) 

2.768 

(0.101) 

0.494 

0.551 

(0.010) 

0.509 

(0.023) 

0.504 

(0.011) 

2.171 

(0.063) 

2.082 (0.129) 

2.105 

(0.083) 

0.715 

0.445 

(0.014) 

0.414 

(0.032) 

0.386 

(0.016) 

1.763 

(0.067) 

1.697 (0.158) 

1.639 

(0.094) 

0.922 

0.397 

(0.021) 

0.331 

(0.049) 

0.299 

(0.021) 

1.627 

(0.090) 

1.129 (0.196) 

1.214 

(0.111) 

1.118 

0.332 

(0.019) 

0.275 

(0.044) 

0.252 

(0.018) 

1.310 

(0.077) 

1.204 (0.176) 

1.028 

(0.085) 

1.303 

0.278 

(0.025) 

0.215 

(0.061) 

0.190 

(0.023) 

1.117 

(0.097) 

0.939 (0.242) 

0.789 

(0.103) 


TABLE XIV. N5 ensemble (a = 0.050 fm, rur: = 429 MeV): Vector form factors at all values for all extraction methods. 


N6 

Ge 

Gm 

0" [GeV^] 

Plat ( 

1.1 fm) 

Summation 

Twc 

I-state 

Plat (1.1 fm) 

Summation 

Two 

I-state 

0.255 

0.681 

(0.011) 

0.644 

(0.012) 

0.588 

(0.012) 

2.618 

(0.096) 

2.797 (0.107) 

2.452 

(0.149) 

0.487 

0.512 

(0.013) 

0.473 

(0.014) 

0.373 

(0.015) 

2.004 

(0.071) 

2.083 (0.083) 

1.717 

(0.124) 

0.701 

0.386 

(0.017) 

0.344 

(0.020) 

0.235 

(0.025) 

1.433 

(0.070) 

1.688 (0.091) 

1.236 

(0.144) 

0.900 

0.363 

(0.026) 

0.331 

(0.031) 

0.173 

(0.036) 

1.341 

(0.102) 

1.754 (0.127) 

1.177 

(0.202) 

1.087 

0.282 

(0.018) 

0.271 

(0.026) 

0.152 

(0.030) 

1.041 

(0.068) 

1.090 (0.101) 

0.724 

(0.158) 

1.264 

0.222 

(0.021) 

0.204 

(0.033) 

0.120 

(0.039) 

0.763 

(0.071) 

0.876 (0.128) 

0.390 

(0.181) 


TABLE XV. N6 ensemble (a = 0.050fm, m-n = 331 MeV): Vector form factors at all values for all extraction methods. 


07 

Ge 

Gm 

Q" [GeV^] 

Plat (1.1 fm) 

Summation 

Twc 

I-state 

Plat (1.1 fm) 

Summation 

Two 

I-state 

0.146 

0.786 

(0.010) 

0.773 

(0.012) 

0.700 

(0.017) 

2.800 

(0.126) 

2.825 (0.139) 

2.886 

(0.259) 

0.282 

0.648 

(0.013) 

0.618 

(0.015) 

0.485 

(0.021) 

2.230 

(0.096) 

2.299 (0.118) 

1.976 

(0.205) 

0.411 

0.554 

(0.015) 

0.521 

(0.019) 

0.357 

(0.026) 

1.929 

(0.086) 

1.964 (0.111) 

1.449 

(0.211) 

0.533 

0.474 

(0.018) 

0.401 

(0.022) 

0.231 

(0.035) 

1.743 

(0.100) 

1.589 (0.111) 

1.059 

(0.239) 

0.650 

0.428 

(0.017) 

0.351 

(0.023) 

0.172 

(0.033) 

1.445 

(0.075) 

1.387 (0.096) 

0.755 

(0.194) 

0.761 

0.388 

(0.019) 

0.327 

(0.029) 

0.176 

(0.033) 

1.249 

(0.070) 

1.266 (0.096) 

0.662 

(0.210) 


TABLE XVI. 07 ensemble (a = 0.050fm, m,r = 261 MeV): Vector form factors at all values for all extraction methods. 
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